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We show that the dc conductance of a quantum wire containing a Luttinger liquid and attached 
to noninteracting leads is given by e 2 /h per spin orientation, regardless of the interactions in the 
wire. It is also shown that weak disorder in the wire results in the temperature- or length-dependent 
corrections to the conductance. The exponents of these dependences are determined by the interaction 
strength in the wire and, in the leading order, are not affected by the presence of the non-interacting 
leads. These results explain recent experiments on quasiballistic GaAs quantum wires. 



I. INTRODUCTION 

For non-interacting electrons, the conductance of narrow ballistic quantum wires connected 
to wide reservoirs is quantized in units of e 2 /h When the effects of interactions are 

included this result is expected to be modified. In particular, when the electrons in the wire 
form a one-dimensional Luttinger liquid (LL) M, the conductance is believed to be Ke 2 /h per 
spin orientation where K is the interaction dependent parameter characterizing the LL. 

For non-interacting electrons K = 1. For repulsive interactions K < 1, and the conductance 
should be reduced. 

A recent experiment on very long GaAs high mobility quantum wires [0 casts doubt on 
this picture. It is known that the same parameter K enters the temperature dependence of the 
impurity correction to the conductance, and using this the authors of @ where able to estimate 
K to be about 0.7 for the electron gas in their wires, implying a conductance reduction of 30% 
in the ballistic limit. The actual reductions observed are only a few percent of e 2 /h however. 

This observation has led us to re-assess the conventional analysis of the charge transport in 
Luttinger wires. In this paper we will argue that the conductance of a quantum wire attached 
to one-dimensional non-interacting leads [which are intended to model the higher- dimensional 
Fermi-liquid (FL) reservoirs] is e 2 /h regardless of the interactions in the wire itselfQ. This is 
because the finite resistance of a ballistic wire is a contact resistance [PHT3"] and comes entirely 



1 The result announced above was obtained by us || simultaneously and independently with San and 
Schulz By now, there are at least three more papers |p~5| — p~T|| which, each in its own way, confirm 
this result. 



from processes that take place outside the wire, where the electrons are not in the LL state (cf. 
Fig. 0a, top). 




FIG. 1. a) Top: a narrow wire containing a Luttinger liquid and connected to Fermi-liquid leads. 
Arrows show a typical electron trajectory giving rise to a finite conductance g = e 2 /h. Bottom: the 
same for a wire with impurities (depicted by asterisks). Dashed arrows show a typical trajectory 
giving rise to a correction to e 2 /h. b) Effective ID model for the system shown in (a) [top]. 

While it has been generally asserted in the literature that g = Ke 2 /h is the correct result for 
interacting electrons (at least for wires longer than the Fermi wavelength), we must remark 
that there have been previous comments supporting the result g = e 2 /h. In their earlier paper 



IS], Kane and Fisher remark that the ac conductance will cross over to the non-interacting 



value at frequencies lower than to = vp/L, where Vp is the Fermi velocity and L is the wire 



length. Matveev and Glazman |J9] make a similar remark when discussing the conductance of 
multi-mode interacting wires. 

Having established the absence of conductance renormalization in pure wires, we will ask 
what would happen if we connect a weakly disordered wire to the leads. In particular, we 
are interested in whether the temperature dependence of the conductance reflects the presence 
of the LL-state in the wire or it is affected by the presence of the leads as well. We will 
show that that the main contribution to the temperature-dependent part of the conductance is 
determined by the local strength of interactions in the wire and is the same as in the absence 



of the leads. This result can be anticipated from the following simple picture | 20fl . Consider 
a weakly disordered wire containing an LL and adiabatically connected to the FL leads (cf. 
Fig. |l|a, bottom). As has been said above, in the absence of disorder, the finite resistance of a 
perfect wire (= h/e 2 ) is entirely due to contact resistance P Hl^l : some of the electrons coming 
from the wide leads are reflected as the channel becomes narrower. This reflection takes place 
outside the wire, where the electrons are in the FL state, therefore, the contact resistance is not 
affected by the interactions in the wire. Weak disorder in the wire gives rise to an additional 
contribution to the resistance. This contribution is determined by the scattering in the wire, 
where electrons are in the LL-state, and therefore this contribution has features typical of an 
LL but not of an FL. 

These two results-the absence of conductance renormalization in perfect wires and the 
temperature dependence of conductance determined by the local strength of interactions in 
disordered wires-resolve the controversy around the experimental observations by Tarucha et 
al. |7|] and suggest that these observations can be considered as an indication of the LL-state 
in GaAs quantum wires. 

This paper is organized as follows. In Sec. II, we prove the absence of conductance renor- 
malization by using the equation of motion for an LL in an external electric field. In Sec. Ill, we 
obtain the same result by using the linear response theory. In Sec. IV, we find the temperature- 



and length-dependent corrections to the conductance due to weak disorder in the wire. Our 
conclusions are given in Sec. V. 



II. LUTTINGER LIQUID UNDER EXTERNAL ELECTRIC FIELD: EQUATIONS 

OF MOTION 

We begin with replacing the original system containing an one-dimensional (ID) wire con- 
nected adiabatically to two-dimensional (2D) leads (cf. Fig. |l|a) by the effective ID system (cf. 
Fig. |l]b). This ID system is an infinite LL which is described by the inhomogeneous effective 
interaction parameter K(x). As is shown in Fig. |I]a, the interaction strength starts from the 
value — 1; e.g., in the left lead, goes through an arbitrary variation in the middle part of 
system ("the wire") and approaches again the asymptotic value of K L = 1 in the right lead. 
The asymptotic value = 1 describes a non-interacting LL and reflects the fact that interac- 
tions in the leads result in the formation of the FL-state, which for the present purposes can 
be viewed as the Fermi-gas state. Similarly, we allow for an arbitrary variation of the density 
wave velocity v (x) in the wire requiring only that it take asymptotic values in the leads, 
which coincide with the Fermi- velocity of the non-interacting electrons. The electrostatic po- 
tential difference applied to the system produces largest electric field in the narrowest part of 
the system, i.e., in the wire. Accordingly, we assume that the electric field E(x) is zero outside 
the wire and undergoes an arbitrary variation in the wire. 

The (real-time) bosonic action for a spinless LL is 

We have chosen to normalize the field so that density of the electrons (minus the background 
density) and the (number) current are given by 

p = d x <f)/^7T, j = -d^jyfk. (2) 

The interaction with an external electromagnetic field is given by (the charge of the electron 
is taken to be — e) 

S mt = ^=fd 2 x {A Q d x 4> - A^} , (3) 

\/7T J 



so the equation of motion for the field (classical or quantum) is 



d t (^d t <j>) - d x (^d x <f) = -L-E(x,t), (4) 



\Kv J \K J y/Trh 

where E = — d x A + d t Ai is the electric field. We assume that the electric field is switched on 
at t = 0, so that E(x,t) = for t < and E(x,t) = E(x) for t > 0. As we see, the problem 
reduces now to that of determining the profile of an infinite elastic string under the external 
force. In this language, (j)(x, t) gives the displacement of the string at point x and at time t, 
while the number current j = —d^fy/ft is proportional to the vertical velocity of the string. 

To develop some intuition into the solution of Eq. (|]), we first solve it in the homogeneous 
case, when K = const, v = const, and E(x) = const for \x\ < and is equal to zero otherwise. 
In this case, the solution of Eq. (01) for late times, i.e., when t > L/v, has the following form 
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for \x\ < L/2; 
for L/2 < \x\ <vt- L/2 
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(5) 



^ for \x\ >vt + L/2, 

where V = EL is the total voltage drop. This solution is depicted in Fig. ^|a. 
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FIG. 2. a) Solution of the wave equation in the homogeneous case for t 
solution in the inhomogeneous case for t S> L/v. 



5L/v. b) Schematic 



The profile of the string consists of two segments (I and II in Fig. |]a) whose widths-equal 
to (vt — L)-grow with time, and of three segments (III, IV and V in Fig. |2]a) whose widths 
are constant in time and equal to L. In segments I and II, the profile of the string <p(x, t) 
is linear in x, and therefore, being the solution of the wave equation, also in £; in segments 
III-V, the profile is parabolic. Outside segments IV and V, the string is not perturbed yet and 
4>(x, t) = 0. As time goes on, the larger and larger part of the profile becomes linear. Thus, 
for late times, the motion of the string can be described as follows: the pulse produced by the 
force spreads outwards with the velocity v, involving the yet unperturbed parts of the string 
in motion; simultaneously, in all but narrow segments in the middle and at the heads of the 
pulse, the string moves upwards with the t- and ^-independent "velocity" dt<f) = KeV/2 ^/nh. 
In terms of the original transport problem, it means that the charge current I — —ej is constant 
outside the wire (but not too close to the edges of the regions of where the electron density is 
not yet perturbed by the electric field) and given by / = Ke 2 V/h. Therefore, the conductance 
(per spin orientation) is g = Ke 2 /h. 

We now turn to the inhomogeneous case. As in the homogeneous case, the profile consists 
of several characteristic segments (cf. Fig. |2|b). In segments III-V, the profile is affected by 
the inhomogeneities in K(x), v(x), and E(x) and depends on the particular choice of the 
x-dependences in all these quantities. In segments I and II however, the profile, being the 
solution of the free wave equation, is again linear in x (and in t). Requiring the slopes of the 
string be equal and opposite in segments I and II (which is consistent with the condition of the 
current conservation), the solution in these regions can be written as 4>(x,t) = A(t — \x\/vj/). 
The constant A can be found by integrating Eq. (^) between two symmetric points ±a chosen 
outside the wire 

-a ^K^) = 7^ J -« dXE(X) = 

t>L, thus A = K L eV/2y/nh. Calculating the current, 



Outside the wire, K(x) = K L and v(x) 
we get g = K^e 2 '/h and, recalling that = 1, we finally arrive at g — e 2 /h. Thus, the 



conductance is not renormalized by the interactions in the wire. The reason for this is that at 
late times (t ^> L/v) the motion of the major part of the string is not affected by the local 
inhomogeneities in the parameters K, v, and L. It is this late-behavior of the solution, which 
one needs to determine the dc conduct anceQ. 



III. KUBO FORMULA: PERFECT WIRE 



We are now going to derive the result that g = e 2 /h regardless of the interactions in the 
wire in a different way, by using the Kubo formula for conductivity. To do so, we adopt a 
simple model, in which the interaction parameter K changes stepwise from the value Kj, = 1 
in the leads (for \x\ > L/2) to the value in the wire (for \x\ < L/2). The electric field is 
assumed to be zero outside the wire but can take an arbitrary variation in the wire. 

The charge current J is related to the field by 

/L/2 r rj/j _ 

dx' / — e-^a^x^E^x'), (7) 
-L/2 J 2lT 

where E w (x) is the temporal Fourier component of the electric field and a^(x, x') is the non- 
local ac conductivity. In the Matsubara representation, cr w (x, x') is expressed via the (imaginary 
time) current-current correlation function by the usual Kubo formula 

a w (x,x') = — / dr(T*j(x,r)j(x',0))e~ 
Tluj Jo 



In the bosonized form, the particle- number current is given by j = —id T <f)/^ and Eq. (|3|) 
reduces to 

a^x.x) = — G^foa^lo-^-e, (9a) 

G%{x, x 1 ) = f driT^x, r)cj>{x\ 0))e~^, (9b) 
Jo 

where G% is the propagator of the boson field (p. The propagator G% (x, x') satisfies the equation 

v(x) 



{ ~ d *(j^ d *) + ^}Gl(x, x') = 5(x - x'), 



(10) 



complemented by the following boundary conditions: i) G%(x,x') is continuous at x = ±L/2 
and x = x', ii) [v (x) / K(x)]d x G% (x, x') is continuous at x = ±L/2 but iii) undergoes a jump of 
unit height at x = x', i.e., 



v(x) 
K(x) 



d x G%(x,x') 



x=x'+0 , . 

= 1. 11 

x=x'-0 y ' 



In addition, we assume that the infinitesimal dissipation is present in the leads, so that 
G°(±oo,x') = 0. (In a real- time formulation this corresponds to outgoing wave boundary 
conditions.) We are interested in the dc limit of (T u (x, x') for which we need to determine the 



2 Safi and Schulz [14| have shown that the ac conductance is renormalized by the interactions in the 
wire in a non-trivial way. 



behavior of G^(x,x') for u — > 0. Note that according to Eq. (^), one needs to know G°(x, x') 
only for \x'\ < L/2. In this region, the low- frequency asymptotic form of the propagator can 
readily be shown to be independent of x, x' and given by G° (x, x') — > iK^/2u ||. Consequently, 
Eq. (Pa| ) gives 

lim <7 w (x, x') = — : — . (12) 

w->0 ft, 

For a static electric field, E w (x) = 2tt5(uj)E(x), and Eq. (0) gives the x- and t-independent 
current 

Ki e 2 r L Kj e 2 

1 ' dx'E(x') = —t—V, (13) 



ft. JO 

from which we see that the conductance is 

9 = K L j. (14) 

Thus the conductance is determined by the value of K L in the leads and does not depend on 
the value of in the wire. Recalling that = 1, we get g = e 2 /h. We have thus confirmed 
the result obtained in Sec. II. 

Notice that in this calculation we have implicitly taken the limit u —>■ while keeping q 



finite [pl| -p3|. The traditional order of limits is q — » before to — > 0. The latter yields the 
Drude formula, which has a divergent dc limit for perfect systems. The former produces a finite 
(two terminal Landauer-Biittiker) dc conductance even for a perfect system. It furthermore 
corresponds to the experimental situation where a static field is applied over a finite region. 



IV. KUBO FORMULA: WEAKLY DISORDERED WIRE 

We now consider the case when a weak disorder potential V(x) is present in the wire. The 
action now acquires the term describing the backscaterring by disorder 

2 r rP 

Si = - dx drV(x) cos(2A; F x + 2^0), (15) 
a J Jo 

in which a is the microscopic length cut-off. In what follows, we will determine the leading 
disorder-induced corrections to directly from the Kubo formula (H) via perturbation theory 
in Si f\. The first non-vanishing correction to the ensemble-averaged propagator is given by 

SGjX, X 1 ) = 4 / dX 1 dX 2 V(x 1 )V(x 2 ) cos(2A; F (x 1 - x 2 ) 
a 2 J 

{(0(X)0(X')Q(X 1 , X 2 )) Q - (0(X)0(X / )) o (g(X 1 , X 2 )) } (16) 

where X = {x,t}, ... stands for ensemble averaging over disorder realizations, (. . .)o for 
thermodynamic averaging and 

Q(X U X 2 ) = e *VWriW(*a)] (17) 



3 The assumption of weak disorder is adequate for the experimental situation of Ref. Q, where L 
was at least ~ 6 times smaller than the elastic mean free path of the unbounded 2DEG and the total 
observed change in the conductance was 1 — 5% of e 2 /h depending on the wire length. 



At this stage, for the sake of simplicity, we choose V(x) in the form of white-noise: V(x) = 
and V(x\)V(x2) = riiU 2 5(xi — X2), where n; is the concentration of "impurities", and u is the 
"impurity" strength. The effective elastic mean free path I (in the absence of the interactions) 
can then be defined as l/l = n\u 2 /a 2 Up, where u>f is the (non-universal) energy ultraviolet 
cut-off (of the order of the Fermi energy)^. We then obtain the correction to the non-local 
conductivity 



TriuJ J-L/2 

[F (x) - F (x)] 



where F^{x) is the r- Fourier transform of the (inhomogeneous) 2kp density-density correlation 
function 



cxp 



F(x,r) = (Q(xr, x0)) = 
-!E,(l-e-^)G°(x,4 (19) 



To evaluate the function F in Eq. fllTf), one needs to know G° {x, x') only for —L/2 < x = x' < 
L/2. Straightforward, albeit lengthy, algebra leads to 

. _ Kw K 2 _e~ L l L * + k + k_ cosh(2x/L^) 

q[X ' X) ~ 2\Q\ \Q\ e^nl - e-V 1 "** ' 1 ' 

where L Q = vw/\u\, t>w is the density- wave velocity in the wire and k± = 1/ifw =t 
We now consider separately the cases of "high" < T < u;f) and "low" (T %/L) 

temperatures. (The quotations marks here are implied to mean that, say, "low" -T case can be 
alternatively viewed as the "long" -L case and vice versa.) 

At "high" temperatures, the second term in Eq. (|20D is exponentially small [oc exp(— [L — 
2x)/Lq] unless u) = 0. The term uj = gives zero contribution to the sum in Eq. (|19D, as can be 



seen by performing the infrared regularization of the propagator (\oo\ — > y |c<}| 2 + m 2 ) and then 
letting m — > 0. Thus, in this case only the first term in Eq. ( ^0|) has to be taken into account. 
This term is precisely the same, however, as in the case of a homogeneous Luttinger liquid with 
the parameter Kyj- Already at this stage the result that the T-scaling of the conductance is 
determined by can be anticipated. The function F is now x-independent and is given by 

27T / y F /3 1 2Xw ^ 



sin ttt / j3 

After analytic continuation u — > iuo — e, the limit c<j — > can be taken. As we saw in 
Sec. Ill, the remaining two propagators in Eq. ( |18|) are x- and x' independent in this limit: 
lim^o G° (x, x') = iK L /2u. The rest of the calculations proceeds exactly as in the homoge- 
neous case. Recalling the result for the conductance of a perfect wire Eq. (|14D , the final result 
for the conductance can be written as 

9 = K h e -h^K h h^] 2{l - K ^\ (22) 



4 Although in the real GaAs system the impurity potential is long-ranged, the expectation is that 
this simplification is not going to affect the final results significantly, as soon as £ is replaced by the 
correct mean free path for a more realistic potential. 



where 



C = S-\/tt sin(7rA'w) 



rn - 



ra + K 



(23) 



Note that A*l enters only the prefactors, while the exponent of T-scaling is determined by A'wf] 
Tracing back the calculations, we can now see the reason for this. The A~l- dependence comes 
from the propagators in the prefactor of Eq. (ET^) which depend on the frequency of the external 
field u. In the limit uo — > 0, these propagators become long-ranged and "know" only about 
but not A'w- The T-scaling comes the propagator entering the function F. This propagator 
does not depend on u. In the " high" temperature limit, it becomes short-ranged and "knows" 
only about A'w but not Al- 

In the case of "low" temperatures, the analysis is much more cumbersome. Referring the 
reader for details to Ref. [20], we give only the final expression the real part of the conductance 



Reg 



— - —2tt 1 (A w )- 

h h £ 1 v w 



2(1-K W ) 



(24) 



Comparing this result with the analogous result for the homogeneous case HU, we see that the 
exponent of the length-dependence is the same as if the leads would be absent^. The L-scaling 
is still symmetric to T-scaling upon replacement T v/L (apart from the additional factor of 
L entering in the combination with the mean free path). Contrary to the homogeneous case 
however, this symmetry exists only if the leads are not interacting, i.e., Kl = 1. 



V. DISCUSSION AND CONCLUSIONS 

In Sees. II and III, we have shown that the conductance of a perfect wire containing a 
Luttinger liquid and connected to noninteracting leads is not renormalized by the interactions 
in the wire. It remains at the noninteracting value g = e 2 /h per spin orientation. In Sec. IV, 
we have shown that weak disorder in the wire leads to the temperature- or length-dependent 
corrections to the conductance, the exponent of which being determined (in the leading order) 
by the interaction strength in the wire and independent of the presence of the leads. 

These results are consistent with the recent experimental observations by Tarucha et al. 
[0, who observed the anomalous temperature dependence of the conductance at lower tem- 
peratures with exponent corresponding to K p ^0.7 but no renormalization of conductance at 
higher temperatures. According to the previous theory, this value of K p extracted from the 
temperature dependence should have implied a 30% renormalization of g from the value of e 2 /h 
at higher temperatures, which clearly contradicts to the data. 

It is certainly tempting to ask whether the observed temperature dependence |7[] really 
comes from the Luttinger-liquid-like behavior of the electrons in the wire or there exists some 



5 Safi and Schulz [25 have found corrections to the leading T-dependent term in Eq. (^) which 
arise from additional scattering at the boundaries between disordered and perfect regions. These 
corrections become dominant only if the interaction is strongly attractive, which is not the case for 
the semiconductor heterostructures. 

6 For electrons with spin, the exponent 2 — 2A"w is replaced by 2 — A~ w — A" w or, if the SU(2) 
symmetry of the underlying Hubbard model is preserved, i.e., iT w = 1, by 1 — A~ w , where K^ 7 are 
the parameters of the charge (spin) parts of the Luttinger-liquid in the wire. 



other plausible explanation for this. One of the possible reasons for the T-dependence might be 
the phase-breaking scattering of weakly-localized electrons at phonons. The counter-argument 
to this assertion is that it would be difficult to explain the anomalous T-dependence of the 
conductance observed in Ref. M (—5g oc T -0 3 ) by the scattering of non-interacting electrons 
at three-dimensional phonons. 

Even within the Luttinger-liquid model, the impurity scattering is not the only possible 
mechanism of the T-dependence. The other one is the scattering of electrons from the non- 
adiabatic openings of the wire into the leadsQ. One can imitate this mechanism in the ID 
model by putting the point scatterers at the ends of the wire, the S'-matrixes of these scatterers 
being equal to that of the non-adiabatic leads. The conductance of such a system is known 
to be T-dependent |p],p6[. However, the T-dependent correction to the conductance is length- 
independent in this mechanism, whereas the observed deviation of the conductance from e 2 /h 
[|7| is definitely more pronounced in longer wires. This gives more support to the disorder 
mechanism considered in Sec IV, although the coexistence of two mechanism is quite plausible. 
We hope to see more experiments in the nearest future providing more checks for the current 
and new theories. 

Finally, we would like to mention that the result of the absence of conductance renormaliza- 
tion by the interactions in a Luttinger-liquid wire discussed in this paper should not be applied 
to chiral Luttinger liquids [27,28] formed at the edges of fractional quantum Hall effect (FQHE) 
systems. The two-terminal conductance of a FQHE system is Qfqhe = ve 2 /h, where v is the 
filling factor. Despite the formal similarity of this expression with the conductance of a homo- 
geneous non-chiral Luttinger liquid (= Ke 2 /h), the physical meanings of these two quantities 
are quite different. In FQHE system, Qfqhe is the Hall (non-dissipative) conductance p9| , p8 



whereas the conductance of non-chiral Luttinger liquids is the longitudinal (dissipative) one. 
Moreover, the finiteness of Qfqhe is not due to the presence of reservoirs (as in the non-chiral 
case). The "renormalization" of Qfqhe from the value of e 2 /h is due to the interactions in the 
bulk, and the bosonic description of a chiral Luttinger liquid |^] is constructed in such a way 
that the correct result for Qfqhe is obtained. 
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7 We thank A. Gogolin for attracting our attention to this issue. 
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CONDUCTANCE DANS UN LIQUIDE DE LUTTINGER 
CONNECTE A DES RESERVOIRS 



Dmitrii L. Maslov et Michael Stone 



Nous montrons que la conductance continue dans un fil quantique qui contient un liquide de Luttinger 
et qui est connecte a des reservoirs non-interactifs est donnee par e 2 /h pour chaque orientation de spin. 
Ce resultat est independant de l'interaction dans le fil. II est aussi montre que la presence d'un faible 
desordre dans le fil resulte en des corrections de la conductance dependant de la temperature et de la 
longueur du fil. Les exposants de ces corrections sont determines par la force de l'interaction dans le 
fil, et, au premier ordre, ne dependent pas de la presence des reservoirs non-interactifs. Ces resultats 
expliquent des experiences recentes sur des fils quantique de GaAs dans le regime quasi-ballistique. 



